Cold bosons in the Landauer setup 
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We consider one dimensional potential trap that connects two reservoirs containing cold Bose 
atoms. The thermal current and single-particle bosonic Green functions are calculated under non- 
equilibrium conditions. The bosonic statistics leads to Luttinger liquid state with non-linear spec- 
trum of collective modes. This results in suppression of thermal current at low temperatures and 
affects the single-particle Green functions. 
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I. INTRODUCTION 

Systems of ultracold Bose gases have recently attracted 
a great deal of experimentalii"— and theoretical attention, 
see Refs. for reviews. A high control over experimen- 
tal conditions, including geometry, density, and interac- 
tion strength, as well as absence of uncontrolled disorder 
allow one to explore new aspects of many-body physics. 
Among experimental achievements, the coherence of non- 
equilibrium Bose condensate was studied based on the 
interference measuremeniji, correlations of density fluc- 
tuations were measured in Refl^, and the distribution of 
the bosons over momenta was explored in Ref. 

Unlike the classic example of a bulk Bose fluid (^He), 
atomic gases are usually realized in optical traps or in 
atom chips where magnetic and electric fields confine the 
system to a geometry with strong asymmetry with re- 
spect to three-dimensional rotations. In many exper- 
imental situations, one deals with arrays of quasi-lD 
systems^iSiiS. These geometrical restrictions strongly in- 
fluence the dynamics as they lower the effective dimen- 
sionality of the system. Indeed, in three dimensions a 
Bose system undergoes a famous Bosc-Einstcin conden- 
sation, its thermodynamic properties are well accounted 
by the mean free theoryii, while its hydrodynamics is 
governed by the Gross-Pitaevskii equation. On the other 
hand, in two dimensions and especially in one dimen- 
sion fluctuations of the order parameter destroy the long 
range order, necessitating a more microscopic treatment. 
In this work we focus on the impact these effects have on 
transport properties of one-dimensional bosonic systems. 

As is well known, a clean one-dimensional (ID) sys- 
tem forms a strongly correlated ground state, so-called 
Luttinger liquid (LL). Though this description holds for 
both, fermionic and bosonic systems, the bosonic char- 
acter begets new properties in Luttinger liquid state. To 
explore these features, we consider Landauer type setup 
shown in Fig[Tl where bosons are trapped in the system 
that consists of two reservoirs connected by a one dimen- 
sional "wire". 






FIG. 1: Two reservoirs with cold atoms are connected by one- 
dimensional trap. The temperatures and chemical potentials 
in the reservoirs are assumed to be different (upper panel); 
interaction strength g{x) as a function of coordinate (lower 
panel) . 



Far-from-equilibrium realizations of Landauer setup 
have been recently studied in the framework of corre- 
lated ID electronic (or, more generally, fermionic) sys- 
tems. The tunneling spectroscopy of the non-equilibrium 
carbon nanotube have been measured in Ref. [T2I . The 
thermal current in the edge state of QHE was studied 
in Ref. [isl . On the theory side, one can distinguish 
several types of non-equilibrium setups, ranging from 
partiall y "'^^'^^ to fully non-equilibrium situationtjiS.. In the 
partial-equilibrium setup, electrons coming from different 
reservoirs have different values of chemical potential and 
temperature. In the case of full non-equilibrium, elec- 
trons in the reservoirs are characterized by an arbitrary 
single-particle density matrix. Remarkably, correlation 
functions of the interacting many-body problem can be 
calculated exactly even in the latter case, and can be cast 
in terms of Fredholm determinants. 

In this paper, we address analogous questions in the 
context of bosonic system. Though at this moment we 
are not aware of any direct experimental realization of 
Landauer-type setup for bosons, the idea seems cxperi- 
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mentally feasible. Indeed, the confinement of bosons to 
ID optical wire has been accomplished in RefH. Since 
the case of partial equilibrium seems to be more natu- 
ral from the point of view of experimental realization, 
we focus on it in the current work. We assume that the 
interaction between the atoms is of a hard core type. 
Inside the reservoirs it plays little role, and we approxi- 
mate the atoms there as an ideal Bose gas. Inside the ID 
"wire" connecting the reservoirs the hard core repulsion 
can not be neglected. We thus describe the system by 
LL with spatially varying interaction parameter g{x), see 
FiglH One of the key features distinguishing this bosonic 
setup from its fermionic counterpart is the absence of 
Fermi surface. In other words, the excitation spectrum 
of particles in the non-interacting regions of the system is 
quadratic. Another interesting realization of excitations 
with quadratic spectrum are transverse spin waves in a 
ferromagnetic Bose gasii. 

The non-linearity of the excitation spectrum is known 
to have a significant impact on the properties of the LL. 
In the fermionic case it leads to a number of interest- 
ing and important effects, see Refs. [I8I - I20I for reviews. 
However, for most characteristics of low-temperature dy- 
namics of correlated electronic systems such non-linearity 
can be discarded. This is done indeed in the case of the 
standard LL model. In the present situation, the spec- 
trum for weak interaction is not linear even in the leading 
approximation. This should be accounted for in the cor- 
responding theory and may be expected to profoundly 
affect the results of our analysis. 

To deal with a non-linear dispersion of the spec- 
trum in the bosonic case, we use the so-called harmonic 
approximation^!. Remarkably, this approach accurately 
describes the system in both the "quasi-condensate" 
(weakly interacting Bose gas) and the LL (sufficiently 
strong repulsive interaction) regimes. Within this frame- 
work, we study the thermal current and single-particle 
Green functions. The latter contain information about 
the density of states and the distribution function of 
bosons, as well as about information about phase- 
coherence correlations that are probed in interference 
experiments^. 



II. BOSONIZATION OF BOSONS 

We begin with the Hamiltonian 

H = Ho + Hi^t, (1) 
that consists of the free part (we set h = 1), 

Ho = -^Jdxi'^dl^, (2) 

and the interaction, 

Hint^ I dxdx'p{x)V{x,x')p{x') . (3) 



Here we define the density of the bosonic field p{x) = 
'i'\x)'ii{x); the bosonic field vj/ satisfies canonical com- 
mutation relations, 



['^ix),¥{x')]^S{x-x' 



(4) 



To analyze the problem we use the hydrodynamic 
approach^!, similar to the bosonization for fermionic sys- 
tems. The term "bosonization" in the present context is, 
perhaps, not optimal since the original system is bosonic 
to begin with. What actually happens is a transforma- 
tion from the original bosonic fields ^^b, ^'^ to new col- 
lective degrees of freedom described by bosonic fields cj), 9. 
The original field operator is expressed in term of the new 
fields a922 23 



i0{x) 



(5) 



Here the field (t){x) is related to the smeared density 
p{;x) = p + Tl{x) (where p is the average density) via 
II(x) — ~dx(t){x)/TT. The collective bosonic fields {(f>,9) 
satisfy the commutation relation 



(x),^(x')] = ysgn(x-x') . 



(6) 



The substitution of Eq. into the Hamiltonian ([T]) leads 
to a hydrodynamic description of ID bosonsii. This is 
a non-linear theory that can be considerably simplified 
by using the harmonic approximation. Expanding the 
Hamiltonian in the fields 9, (f>, one keeps terms only up 
to the quadratic level, which yields 



2m 



dx 



^{d,p)^ + p{djf 



(7) 



The harmonic approximation allows us to account for the 
non-linear (in the present case quadratic) spectrum of the 
low energy sound mode in the collective theory. Its valid- 
ity is restricted to the low-energy (large-density) regime, 
T <C maxjAjpg}. Here A = p^/m is the bosonic coun- 
terpart of Fermi energy, and g is the interaction strength 
introduced below. 

Clearly, the harmonic approximation is not an exact 
theory. By neglecting the interaction between low-energy 
modes (represented by terms of higher order in p and 9 in 
the Hamiltonian), one discards relaxation precesscs that 
are important, in particular, for thermal equilibration, 
drag, and thermoelectric effeciii^i^'^ir— . We will assume 
that our "wire" is not too long, so that neglecting these 
processes is justified. In that situation the harmonic ap- 
proximation is sufficient to describe thermal transport 
and tunneling spectroscopy in the system. 

In terms of the collective bosonic fields 9 and (/) we 
obtain 



Ho = 



dx 



1 



dx 



8TT^mp{x) 



2m \ 



(8) 
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Here we model the interaction with a short range poten- 
tial, 



V{x, x) = g{x)5{x — x') . 



(9) 



Let us note that Eq. (jH)) in fact corresponds to the micro- 
scopic model (jni) in the limit of weak interaction only. In 
the local interaction model the large-g limit yields 
the Tonks- Girardeau gaa^^ that can be mapped onto 
a free-fermion model (characterized by the LL param- 
eter K = 1). However, it makes sense to consider g 
in Eq. ([8]) as a phenomenological parameter of underly- 
ing LL model, which allows us to go beyond the Tonks- 
Girardeau limit. On a microscopic level, this corresponds 
to the replacement of the delta-like repulsion (|9]) by a 
finite-range hard-core interaction. 

Equation ([5]) is the Hamiltonian of a ID LL for inter- 
acting bosonie system in the harmonic approximation^. 
As one clearly observes, the Hamiltonian contains 
a fourth-order spatial derivative [the term (9^0)^], at 
odds with the standard fermionic LL model that con- 
tains only a second spatial derivative [{dx(j)Y]- For 
large values of the interaction constant, pg » p^/2m, 
the fourth-derivative term is relatively small, leaving us 
with the standard LL with the interaction parameter 
= 'ir^p/2mg. When the interaction is weakened, 
pg ^ p^/2m, the standard (linear-spectrum) LL descrip- 
tion is valid for the lowest temperatures only, T pg. 
The LL parameter becomes larger for weaker interaction 
and tends to infinity in the limit of free bosons but the 
region of applicability of the standard LL theory vanishes 
in this limit. 

In the general case, one should use the full theory 
(rather than the standard LL theory). The correspond- 
ing spectrum of bosonie excitations uj{q) is non- linear— 
, interpolating between quadratic (ui = /2m for non- 
interacting bosons) and linear (w = uq for strongly in- 
teracting bosons; u being the sound velocity). Finally, 
we mention that the Hamiltonian (jH)) depends on the 
mean bosonie density, p{x). To find the profile of the 
bosonie density, one needs to go beyond the LL descrip- 
tion and solve the non-linear hydrodynamie equations, 
see Appendix IB 41 

To deal with the non-equilibrium conditions, we use 
the Keldysh formalism. The fields on the upper and lower 
branches are labeled by -I- and — respectively. It is con- 
venient to perform a rotation in Keldysh space 



0,0= (0+±0_)/x/2, (10) 

e,e = {e+±e^)lV2. (ii) 

where we will refer to (0, 6) as classical and to (0, 0) as 
quantum components^. We then find that the system is 
described by the following action 
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(12) 



where we have defined the vector $ = (0, 9, 4>, 9) . 
The inverse propagator has a standard structure in the 



Keldysh space 







(13) 



where each component is a matrix in the 0, space. The 
inverse retarded propagator is given by 



2-K 2m"a:/^(-'a 



(14) 



with 



- ^ + d^-d, ~ 8mg 



p- 



p- 



P 



and a;-|_ =0^-1-10. The advance component is related to 
the retarded one by complex conjugation, {D~-^)°-{oj) = 
[{D-^Y{uj)Y. The Keldysh component (D"!)^ carries 
information about the distribution functions in two reser- 
voirs: right-moving modes coming from the left reservoir 
have the temperature Tr, while the left-moving modes 
coming from the right reservoir are characterized by the 
temperature T^. 

To simplify the analysis from now on we consider the 
large density limit, approximating the mean bosonie den- 
sity p{x) by a constant. Variation of the action with re- 
spect to the classical components of the fields 9 and 
yields the saddle point equations 
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^dj^ + ^dl9^ = 



27r 



0. 



27r 



2ra 



(15) 



These two equations can be conveniently combined into 
the wave equation 



2pg{x) 



1 



4™^ 



J{uj, x) 



0, 



(16) 



where we have introduced J ~ Tr^^dx9. Equation ([2 
describes propagation of plasmons inside the wire with 
dispersion that varies in space (as a result of the vari- 
ation of g{x)). Due to these variation, a plasmon may 
experience scattering but the total number of plasmons 
is conserved. To make this conservation explicit, we mul- 
tiply Eq. ([TS]) by J* from the left and subtract a complex- 
conjugated equation. This yields the plasmons conserva- 
tion law 



dtQ + dxJ = 



(17) 



which has a form of the continuity equation that states 
that the charge 



Q = JdtJ* - J*dtJ 
is carried by the current 



(18) 
(19) 
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FIG. 2: The boundary between the regions with different in- 
teraction constants. Propagating and decaying waves on both 
sides of the boundary are shown. 



FIG. 3: The scattering of plasmons at the boundary between 
regions with different values of interaction can be described 
by scattering matrix S. 



In a region with a constant interaction strength g 
Eq. yields the Bogolubov's excitation spectrum for 
the acoustic phonons, 



2P5 



4m^ 



(20) 



This dispersion relation has four solutions, resulting in 
oscillating and exponentially decaying (growing) waves 
(see Fig.E]), 



where 



(21) 



= V2m^J-g+ ^yg^ +ujym^, (22) 
= V2m^Jg+ ^g^ +oj'^/m^. (23) 



We consider now the situation when the interaction 
strength changes from one value to another in a boundary 
region, see Fig.O If we consider the solution not too close 
to the boundary, the exponentially decaying components 
can be neglected, and the propagating waves are related 
via a scattering matrix. To take into account the different 
velocities of propagation, we define a ~ ^/uA, b — y^B, 



the transmission and reflection amplitude of such a bar- 
rier, we derive matching conditions for the amplitudes at 
the boundary. For this purpose, we integrate Eq. (jl6p 
over a small region around the boundary. This leads to 
a requirement that J, J', J", and 2pgJ' — J'" /Am are 
continuous at the boundary [which implies that J'" has 
a jump equal to 8mpJ'{g2 — 5i)]- These four conditions 
allow us to find the amplitudes A2, Bi of the outgoing 
waves and Ci, D2 of the decaying waves for given am- 
plitudes Ai and B2 of the incoming waves, and thus to 
establish the scattering matrix. We obtain the transmis- 
sion 



t 



2 / ((?i +ipi)((?2 
'9i'?2- 



W2) 



Z y """" {qi - ipi){q2 - ip2) 
X [8TT'mpp2{gi - 52) - {pi + P2){ql + pI)] , 

and the reflection amplitude 
1 52 + iP2 



(26) 



Z q2- ip2 
+8Trmp{piqi +P2'?2)(ff2 - 5i) 



- [pi +P2){qi + ip2){q2 - ipi){qi - 92) 

(27) 



^ is a sound velocity in corresponding ^^^^^ j^^^^ introduced the notation 



where Uq 

region. It is easy to verify that coefficients a and h in 
different regions (see Fig. [3]) are related 



'02 \ Ml 
I \b2 

through a unitary scattering matrix, 



(24) 



(25) 



Z = {pi+ P2){qi + ip2){q2 + ipi){qi + q2) 

+8Tnnp(j)iqi - P2q2){g2 - 9i)- (28) 

For our model with non-interacting reservoirs, we now 
consider the case of the scattering between interacting 
and non interacting regions, i.e. gi = and g2 = g. In 
this case the transmission amplitude is a function of a di- 
mensionless parameter s = uj /lupg with the asymptotics 



with \t\ = \t'\ and |r| = |r'|. 

The transmission coefficient has to be calculated for 
a particular realization g{x) of the interaction in the 
boundary region. In the limiting case of an adiabatic bar- 
rier, when the interaction changes smoothly on the scale 
of the plasmon wave length, one finds the ideal trans- 
mission, \t\ = 1 and r ~ Q. We focus on the opposite 
case of a sharp-step barrier, with interaction constant gi 
to the left of the boundary and (72 to the right. To find 



_23/4si/4^ for s < 1 
1, for s > 1. 



(29) 



III. 



KINETICS OF ID BOSE FLUID: THERMAL 
CURRENT 



In the preceding section, we have found the plas- 
mon transmission coefficients at the boundaries between 
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FIG. 4: Distribution functions for transmitted and reflected 
plasmons inside tlie wire 




FIG. 5: Total plasmon transmission for a system of two iden- 
tical sharp barriers as a function of s = u}/2npg. 



the interacting region and the reservoirs. Supplement- 
ing these results with the boundary conditions on the 
distribution functions of plasmons in the reservoirs, we 
straightforwardly calculate the distribution functions of 
the right and the left moving modes {Bj^/j^) inside the 
wire, see Fig. |31 

For a wire longer than the thermal wave length of plas- 
mons, the Fabry-Perot-type plasmon interference can be 
neglected, and one finds 



Br- 



Ti 



7^l7^2 



-B 



(0) 



Br 

7^27^l ^ 



(0) 



72Jll_ 

-7^l7^2 
rl7^2 



-B 



(0) 



Br 

7^27^l « 



(0) 



(30) 



where = coth(a;/2Tfl), B'^'' = coth(a;/2TL), Tl and 
Tfl are temperatures of the left and the right moving col- 
lective modes, 71 = T2 = 1^2 P are the transmission 
coefficients of the left and right barriers, and TZj = 1 — Tj 
are the corresponding reflection coefficients. . 

Using the plasmon distribution function one can cal- 
culate, in particular, the thermal current 



lE = ^J^^duju[B''i:{u)-BZ{u)]. 
Substituting Eq.dSD]) into Eq.dSJ), we find 

dojij[Bl{ij) - BU^)]Ttot{uj) . 



Here Ttot is the total transmission coefficient 



Ti(a;)7^(a;) 
l-7^l(w)7^2(a;) 



(31) 



(32) 



(33) 



it is shown in Fig. [5] for the case of two sharp barriers 
characterized by the transmission amplitude p6|) . 

Substituting Ttot for the case of two sharp barriers into 
Eq. (j32p and performing the frequency integration, we 
find 



where 



Ie=J'e{Tl)-J'e{Tr), 



, for T > 271 pg 



(34) 



(35) 



In the limit of small temperature difference AT — ~ 
Tr wc have 

'" = 8^i '" sinh^l/2T '^"^"^' ^''^ 

where T ~ 2i±Ik Eg. ([55]) can be cast in the scaling 
form 



(37) 



The function / entering Eq. ((37)) can be calculated nu- 
merically and is plotted in Fig. [HI It has the following 
asymptotic limits: 



/(^) 



as^/^, for s < 1 , 
I, for s>l. 



(38) 



where a ~ 1.003. 

At relatively high temperatures, T 3> pg, we repro- 
duce the thermal current of non-interacting bosona^, 
see Appendix This result can be considered as a 
thermal-current counterpart of the Landauer quantiza- 
tion of charge conductance: the numerical coefficient in 
the second line of Eqs. ([55)) and (|55)) is fully univer- 
sal in ID geometry and does not depend on the form 
of the spectrum, interaction strength, carrier statistics 
(fermions vs. bosons), etc. The only condition is the 
absence of back-scattering of plasmons. In the low- 
temperature regime, wc find that the thermal current 
is suppressed due to reflection of the bosons on the 
boundary between interacting and non-interacting re- 
gions. This r^/^ suppression of the thermal conductance 
distinguishes the bosonic setup from its fermionic coun- 
terpart. 



IV. PHASE COHERENT DYNAMICS: GREEN 
FUNCTIONS 

We now proceed with the analysis of the single-particle 
Green functions of the original bosons. 



G>(x,T) = -*(*B(x,i)«'L(0,0)) 



(39) 
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FIG. 6: Scaling function f{s) governing the dependence of 
the heat current on s = T /2npg, see Eq. (|37|) 



that carry information about spectral properties (den- 
sity of states and distribution functions) of the system. 
The results for the non-interacting case are well known; 
for completeness we present them in Appendix [Xj Our 
approach allows us to analyze the Green functions in 
a broad range of parameters, with the only assump- 
tion being T <^ ma,x{pg, p'^ /2m). To relate it to well- 
accepted terminology in the field^i^, the harmonic ap- 
proximation allows us to describe the system in the "Dc- 
cohcrcnt Quantum" , "Quasi-Condensate" , and "Tonks- 
Girardeau" regimes, and in the strong-interaction LL 
regime, pg 3> T, p'^ /2m (which is not realized in the 
delta- interaction model considered in Ref. [3^|33[ ) , as well 
as in crossovers between them. 

Using the representation of the boson creation and 
annihilation operators in terms of the collective field, 
Eq. (O, we write the Green function as a correlation 
function of the harmonic fluid. 



G^(x,t) 



-z(rKP^/'(x,T)e^«"T(^^^)pi/'(0,0)e-'^'±("'°)). 

(40) 

It is convenient to represent this correlation function 



as 



G^{X, t) 



-ip 



d 

dai 



d 
da2 



92 



X e 



Here we have defined a four-component 

-1 d ' 



daida2 ^ 

(41) 

source" vector 



27^p^/2 dx 
1 

T2 
±1 d 

27rpV2 dx 



4,— CJ.X 



1 



b(x -X)- ^(x 
e~"^^ai5{x -X)- a2(5(x) 
5{x^X) + 5{x) , (42) 



and it is understood that one should set the sources 
ai,2 = after the derivatives in Eq. (|4ip have been eval- 
uated. 



Since the action ([12)) is Gaussian, the functional in- 
tegration over fields 9 and (f> can be easily performed, 
yielding 



. ' . f d d \ 92 
' ^ \dai da2 J daida2 



X2 



Thus, the problem of calculation of the Green functions 
has been reduced to the calculation of the bosonic prop- 
agator 



D 



D" 



(43) 



We now focus on the case of coinciding spatial points 
{X = X') deeply inside the interacting part of the wire 
(for X X' see Appendix IB 31) . Expanding the bosonic 
density operator up to second order in (/), and calculating 
Gaussian functional integral over the bosonic fields (see 
Appendix IB 31 for technical details), we obtain 



G^{t) = -ipil-^fir) 



(44) 



Here we have defined the pre-exponential factor 



$<(r) 



27rp Jo "ipgq + q^/m 



(45) 



I ■ q \ *9 . 

X iio sm LOT cos ujt ± sm lot =F 2lo cos lot 

' 4to / 2to 



and the exponential factor 



2p 



d^^-i^o/A 9 + q'^/8rnp 
27r q{g + q'^/4mp) 



{B^ +B^){1- cos ujt) ± 2i sine 



(46) 



In these formulas q should be understood as related to uj 
via the dispersion law (P^ . 

It is instructive to compare our results for the Green 
functions of a bosonic system in a partial non-equilibrium 
state, with their fermionic counterparts^^. It is seen that 
the bosonic results differ in two respects: (i) the appear- 
ance of a pre-exponential factor and (ii) a more compli- 
cated form of the exponential factor reflecting the non- 
linear character of the spectrum of collective excitations. 

In fact, Eq. ([^^ interpolates between the standard- 
LL results (applicable in the limit of strongly interacting 
bosons) and Ea. (|Bl[) which is valid for free bosons. The 
characteristic energy scale for the crossover between these 
two regimes is set by the interaction (ujq = pg). For en- 
ergies well below loq, the model behaves as a standard 
LL system (with a linear plasmonic spectrum). In this 
case the pre-exponential factor (E>i is small (of the order 
of T/loq ^ 1), and therefore can be neglected. The ex- 
ponential factor $2 in this limit turns into the LL result 
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of Rcf. [Tsl (which, of course, reduces to the conventional 
LL formula in the equilibrium case). In the high-energy 
limit (max{r~^,T} ^ ujo) the spectrum of excitations 
is quadratic, and one recovers the free-boson results, see 
Appendix |^ 



V. SUMMARY 

In this work, we have studied a system of bosonic 
atoms in a Landauer setup, subject to temperature and 
chemical potential difference. We have developed a non- 
equilibrium bosonization approach that describes the sys- 
tem within the harmonic approximation. This approach 
ignores the interaction between different collective modes 
but takes into account the non-linear dispersion of their 
spectrum. 

We have studied the plasmon propagation in a two- 
terminal setup formed by two non-interacting reservoirs 
connected by an interacting ID "wire" . The plasmon 
back-scattering is controlled by the dimensionless param- 
eter s ~ i^/pg, where w is the plasmon frequency (whose 
characteristic value is set by the respective tempera- 
tures of the reservoirs), p is the density of bosons, and 
g the interaction strength. At large s, back-scattering 
is suppressed, and the thermal current acquires a uni- 
versal value (that depends neither on the spectrum of 
the original particles nor on their statistics) . In the low- 
temperature regime, reflection of plasmons is strong, and 
the thermal current exhibits a T^/"^ suppression com- 
pared to the universal value. As another application of 
this formalism, we have calculated the bosonic single par- 
ticle Green function in the same non-equilibrium setup. 
The result interpolates between the conventional (linear 
spectrum) Luttinger-liquid and free-boson limits. 

The approach developed in this work can be used to 
study other properties of the system, e.g., higher cor- 
relation functions or non-steady-state characteristics. A 
more ambitious perspective will be development of a gen- 
eral non-equilibrium bosonizcd theory of non-linear LL 
(formed by interacting fcrmionic or bosonic particles) in- 
cluding both non-linear spectral dispersion of plasmon 
modes and their coupling. 



VI. 



ACKNOWLEDGMENTS 



We thank N. Davidson, L. Khaykovich, and J. 
Schmiedmayer for useful discussions. Financial support 
by German-Israeli Foundation, the Israel Science Founda- 
tion, and the Minerva Center is gratefully acknowledged. 



Appendix A: Transport properties of free bosons 

In this Appendix we summarize basic properties of 
non-interacting bosons, using the original description. 



1. Landauer approach: Particle current 

For the case of non-interacting bosons the particle cur- 
rent can be straightforwardly found within the Landauer 
approach: 

poo 

1= / deiy{e)v{e)[NL{e) - Nn{e)] , (Al) 
Jo 

where v is the density of states, v the velocity, and Nj 
the distribution function in the corresponding reservoir. 
Using the relation v{e)v{e) = l/27r, one finds the relation 



/ = 



(A2) 



Assuming Bose distributions with the same temperature 
T and two different chemical potentials pL, PR. and per- 
forming the integration, we obtain 



^"^2^ ll-eA^H/T 



(A3) 



We note that the particle current is determined by the oc- 
cupation numbers in both reservoirs at the bottom of the 
band (at e = 0). In the linear response regime, expand- 
ing the current in the difference of chemical potentials, 
we find 



I=^{pL-PB)NBie^O). 
Zn 



(A4) 



To relate the occupation number at the bottom of the 
band to macroscopic parameters of the problem we com- 
pute the particle density 



deiy{e)NB{e) = p. 



For \p\ <^T one finds 



p = T, 



thus leading to Nsie = 0) = /mT and 



T ^ P ( 



pr) ■ 



(A5) 



(A6) 



(A7) 



We note that the particle current of bosons is enhanced 
by a large parameter p^/mT as compared to the particle 
current of fermions subjected to the same difference of 
chemical potentials. The large conductance is reminis- 
cent of superfluidity of Bosc condensate in higher dimen- 
sions. 

On a more formal level, the appearance of the "Fermi 
energy" p'^/m entering the large factor p^/mT indicates 
that the particle current can not be calculated within 
the harmonic approximation. In the interacting case its 
calculation would require the use of a full non-linear hy- 
drodynamic theory, see Appendix IB] below. 
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Thermal current 



2. Non-interacting bosons: Bosonized description 



The thermal current is given by 
de 



Jo 

Performing the integration one finds 



(AS) 



(A9) 



We note that the heat current in the system of non- 
interacting bosons coincides with the heat current of ID 
free fermions^. Moreover, this result is universal and 
does not depend on the shape of single particle spectrum 
of elementary carriers, either fermions or bosons^. This 
universality survives the adiabatic switching of interac- 
tion (again for both fermions and bosons), when back- 
scattering of plasmons is negligible, but is ultimately vi- 
olated in the generic case [see, in particular, Ea. (|37|) ]. 



Appendix B: Green functions, G< 

In this Appendix, we present details of the calculation 
of the bosonic Green function. We begin by consider- 
ing the case of non-interacting bosons, first in the origi- 
nal formulation and then in the within the bosonization 
framework. Then we calculate the Green function in the 
interacting case. 



The Green function within the bosonization framework 
is given by Eq. (j44|) . Substituting the spectrum of free 
bosons in Eq. we get 



> 1 [rn f°° doj 



± — sin UT — cos LOT 
2 



(i sin LOT — — cos ujt) 
(B5) 



Similarly, Eg. (1461) yields for the function in the exponent. 



1 



(1 — cos,ojt)B^ 



±i sincjr 



(B6) 



Let us emphasize that the integral over frequency in 
Eq. (|B5|) converges (as opposed to the logarithmically di- 
vergent integrals in the LL case). The resulting function 
$2 is actually small. Therefore, one can expand the ex- 
ponent in Eq.(j33]). Combining all terms together, one 
finds 



(B7) 



To compare this with the exact result, we take Eq. (|B1|) . 
subtract its value at r = (density), and then consider 
the limit /i — > 0. The result is in full agreement with 
Eqs. 113), and (|B6|) . 



1. Non-interacting bosons: Original description 



3. Interacting bosons, Bosonized description 



We consider the non-interacting ID bosons and cal- 
culate the Green functions for coinciding spatial points, 
X = X' (the case oi X ^ X' treated similarly). Equa- 
tions ^ yield 



G|(r) 



2tt 



2T 



(Bl) 

In the energy representation it reproduces well known 
results 

G<{uj) ^ ~2TTiiy{uj)NB{uj - p), (B2) 
G>{io) ^ -2T:iv{uj)[l + Nb{uj - p)], (B3) 



where Nb is the Bose distribution function and 

e{uj) 



(B4) 



is the density of states of non-interacting bosons. To 
compare with the bosonized theory, we consider the high- 
density (/i — 0) limit. Next, we show how these results 
can be derived using the bosonization approach. 



We now consider the case both coordinates X and X' 
are located inside the interacting part of the system. In 
this case the Green function depends on the x = X — X'. 
For the problem of the sharp barrier one finds 

Here the pre-exponential 



duj 



'ipgq + /m 



B'fj^ iuj sin 



and exponential factor 



duj g + l%mp 



2p Jo 27T q{g + q^/Amp) 
B^{l~cos^R)±ism^R + {R^ L)] , (B9) 



where Cr/l — ^T^qx and / is an average value of the par- 
ticle current fiowing through the system, see Appendix 
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IB 41 It is given by the mean value of (dxd), which is 
determined by a non-linear hydrodynamic equation. 

We note that calculating the correlation function of the 
bosonic fields we neglected the modulation of the mean 
density of bosons. It is possible to generalize the har- 
monic approximation and to allow the modulation of the 
mean density in space. We now present general formulas 
applicable in the case of spatially varying mean density. 

To calculate the correlation functions of bosonic fields, 
one needs to find the inverse of the operator 



D-'^D = 1 



(BIO) 



Employing Eq. ([Ti|) . we find that components of the 
matrix correlation function D satisfy the following equa- 
tions: 



m 

LdxD'gg{uj,x,x') = {27T)^K-\x)d~H, 
where we have defined an operator 



L = -K'-^x) . 

m 



(Bll) 



(B12) 



In order to find the components of the correlation func- 
tion D, we construct the scattering states Xq,Tjix), char- 
acterized by the momentum q and index rj that labels the 
reservoir from which the state was "emitted" . The scat- 
tering state wave function satisfy the following equation 



K-'x'',{x)^io^,x'g{x). 



(B13) 



and similarly for other components. 

4. Mean density current for interacting bosons 

To find the particle current one needs to calculate the 
expectation value of the operator 



-dJix) . 



(B17) 



In the presence of a weak external potential U {x, t) the 
linear response theory predicts 



{I{x,t))^dt j dx'dt'Dl^{x,t-x',t')d^U{x\t'). (B18) 

We now analyze this expression in different limits. If the 
interaction between bosons in the leads is finite (we 
now relax the assumption = we used throughout the 
manuscript), it sets the energy scale pg^- For energies 
below this scale the spectrum of collective modes is linear, 
and one restores the known LL resulli^^^ 



(/) = KrV/h . 



(B19) 



Here Kr ~ \/ ■n'^p/2mgr is the LL parameter in the leads, 
and 1^ is a difference between the chemical potentials in 
the leads. In the absence of interactions in the leads, the 
value of current in this case is unaffected by the interac- 
tion inside the system, as in the fermionic case^"— . 

For the case of free bosons [gr ~ g ~ 0) Eq. ljBlSp 
yields 



(B20) 



In terms of the scattering states, the correlation func- 
tions of the bosonic fields can be written as follows: 



i?;/;(u;,x,x') = ^p'■/^ 



(B14) 



Dl';{u:,x,x') = 2^^u^-^'^-^V^'\ 
Dli\u,x^x') = -{2^fk-\x')d-^d-^V^/\ 
Here we defined 



(B15) 



The Kcldysh component can be constructed from the re- 
tarded and advanced one, by imposing the "partial equi- 
librium" condition in each direction of plasmon propaga- 
tion. 



V x'I,{x)x7(^^')Br,{^)Ku^'-uj^), (B16) 



q>o,v 



implying that the "conductance" diverges in the d.c. 
limit (oj — >■ 0). This divergence signals that the harmonic 
approximation is not a suitable framework to calculate 
the particle current. Indeed, while developing the har- 
monic approximation we took the /x — >■ limit, assuming 
that characteristic energies of the collective excitations 
are much greater than the chemical potential. While 
this assumption is valid for the thermal current and the 
single-particle Green function it docs not hold for the 
particle current of the free bosons. To cure this prob- 
lem, one should restore the low energy cut-off, replacing 
Lo with fi. Doing so and using Eq. (jA6p . one recovers the 
conductance of non- interacting bosons, Ea. (|A7[) . 

While the limits discussed above give us some idea 
about the particle current, this problem remains to be 
solved for bosons that do interact in the wire {g 0), but 
do not interact in the leads {gr = 0). In particular, it re- 
mains to be seen whether the value of particle current in 
this case is affected by interaction inside the system. To 
answer this question, one has to go beyond the harmonic 
approximation and to resort to the full (non-linear) hy- 
drodynamic description. 
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